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Abstract 



The rainbow truncation of the quark Dyson-Schwinger equation is combined 

with the ladder Bethe-Salpeter equation for the dressed quark-photon vertex 

to study the low- momentum behavior of the pion electromagnetic form factor. 

^ ' With model gluon parameters previously fixed by the pion mass and decay 

^^ . constant, the pion charge radius r^ is found to be in excellent agreement with 

the data. When the often-used Ball-Chiu Ansatz is used to construct the 

^D I quark-photon vertex directly from the quark propagator, less than half of r^ 

Q.^ I is generated. The remainder of r^ is seen to be attributable to the presence 

0> ' of the p-pole in the solution of the ladder Bethe-Salpeter equation. 
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I. INTRODUCTION 

For timelike photon momenta Q^ in the vicinity of the p-meson mass-shell, the pion 
charge form factor Ft^{Q'^) will exhibit a resonant peak associated with the propagation of 
intermediate state p-mesons (we ignore a possible small effect due to p — uj mixing). That is 

p^^Q') - ,^.f^;'"^ p, , (1) 

Qp [Q^ + m^ - inip Tp) 

in the Euclidean metric (Q^ < corresponds to the timelike region) that we use throughout 
this work. Here rn?p/ Qp is the p — 7 coupling strength fixed by the p -^ e^ e~ decay, f^p^^r 
is the coupling constant for the p —>■ nir decay, and Tp is the p width which is principally 
due to the latter process. A long-standing issue in hadronic physics is the question of the 
extent to which F-j^^Q"^) at low spacelike Q^ can be described by the p- resonance mechanism. 
This is an essential element of the vector meson dominance (VMD) model which is one of 
the earliest field theory models to be successful in a point coupling description of aspects of 
hadron dynamics. In the form of VMD where p — 7 coupling is described by the contraction 
of the two field strength tensors p^'^F^,^, the pion charge form factor is produced as 



The non-resonant first term arises from the photon coupling to the charge of a point pion. 
The resonant second term arises from p — 7 coupling, implements all of the Q^-dependence, 
and vanishes at Q^ = in accord with gauge invariance. The width Tp is nonzero beyond 
the threshold for tttt production only, and the form factor is real for Q^ > — 4m^. In 
this model the charge radius r^ = —6F^{0) comes entirely from the resonant term and is 
Ggp-KTr/i'm'lgp) = 0.48 fm , which compares favorably with the experimental value 0.44 fm . 
See Ref. [0] for a recent review of the pion charge form factor and VMD in this form. 

In terms of QCD, where the pion is a qq bound state, the content of Eq. (|^) cannot 
provide a realistic picture for F^^^Q^) at (low) spacelike Q^. The photon couples only to 
the distributed quark currents in the pion, the vector meson bound state is not a well- 
defined concept away from the pole, and the question of a resonant p contribution has to be 
addressed within the dressed quark-photon vertex. This is the topic we explore in this work. 
Given that one can find a convenient (and necessarily model-dependent) representation 
of the quark-photon vertex involving direct and resonant parts, the direct coupling will 
necessarily produce a distribution F^*'''((5^) to replace the first term in Eq. (Q) and there 
will be a corresponding contribution to the charge radius. The remaining contribution must 
differ from the second term of Eq. (^ because gpTrn, dp, "^p and Tp are well-defined only for 
the physical on-shell vector meson bound state. A measure of the ambiguities involved is 
provided by a recent study that modeled the underlying qq substructure of p and tt. Under 
the assumption that the p Bethe-Salpeter (BS) amplitude is applicable also at Q^ = 0, it 
was found that gp^^{Q^ = 0) ~ gpTrnf^ H]- On face value, this suggests that only about 50% 
of r^ would be attributable to the p-resonance mechanism. 

In this work we obtain the quark-photon vertex as the solution of the inhomogeneous 
Bethe-Salpeter equation (BSE) in ladder truncation with the dressed quark propagators 



taken as solutions of the Dyson-Schwinger equation (DSE) in rainbow truncation. Such 
a procedure automatically incorporates the pole structure of the vertex corresponding to 
the vector meson spectrum of the homogeneous BSE. A previous exploratory study of 
the coupled DSE-BSE for the quark propagator and the quark-photon vertex employed a 
simple infrared dominant form of the effective gluon propagator in order to utilize closed 
form expressions for the dressed quark propagators and the resulting BSE kernel. That 
study was not carried far enough to draw implications for the pion form factor. 

We employ a more realistic model for the effective quark-antiquark coupling that has 
recently been shown to reproduce the pion and kaon masses and decay constants ^ as well 
as the masses and decay constants for the vector mesons p, (f) and K* to within 10% 0. 
There, as well as here, the quark propagators are consistently dressed through the quark DSE 
using the same effective gluon propagator, which ensures that the vector Ward-Takahashi 
identity is obeyed and electromagnetic current is conserved. The model parameters are all 
fixed in previous work |^ and constrained only by m.„, mx, fn and {qq). 

Next, we use this quark-photon vertex to investigate the low momentum behavior of the 
pion charge form factor. The results are compared to those from the Ball-Chiu (BC) Ansatz 
for the vertex which is commonly used in such studies. In comparison with the results from 
the BSE solution, use of the BC Ansatz (and variations thereof) underestimates the charge 
radius. The BC Ansatz does not have the vector meson pole that is naturally present within 
the BSE solution. We explore the extent to which the BSE solution for the vertex can be 
reproduced by the addition of a suitable p-pole term to the BC Ansatz. We then assess the 
performance of this "resonant-improved" Ansatz for the quark-photon vertex. 

In Sec. II we review the formulation that underlies recent studies of the pion form factor 
within a modeling of QCD through the DSEs, and discuss the need for a dressed quark- 
photon vertex in such models. The BSE for the quark-photon vertex is described in Sec. Ill 
and the vector meson pole contributions are outlined. Also described there are: the model 
used for the effective coupling in ladder truncation, the basis of covariants used, and the 
comparison of the numerical solution with the resonant-improved BC Ansatz. Our analysis 
of the low-momentum behavior of the pion form factor and charge radius is presented in 
Sec. IV, and concluding remarks are given in Sec. V. 

II. THE PION ELECTROMAGNETIC FORM FACTOR 

A number of works have treated the electromagnetic form factors of pions [§|,0, 
kaons @,0, and more recently vector mesons [^ and nucleons [|ll|] in the course of QCD 
modeling via truncations of the DSEs. Using dressed quark propagators and bound state BS 
amplitudes, form factors can be calculated in impulse approximation, as depicted in Fig. |l| 
for the 77r7r vertex. Such work proceeds most easily in Euchdean metric {7/i,7i/} = 25^,^, 
7^ = 7^ and a ■ b = J2i=i ^i^i- ^^ ^^^ space of color, flavor and Dirac spin, we denote by 
r^(p;Q) the quark-photon vertex describing the coupling of a photon with momentum Q 
to a quark with initial and final momenta p — Q/2 and p + Q/2 respectively. With this 
notation, the 'Jtit: vertex takes the form 



= -2 1^ Tr [V^ik^; -P+) 5(g+_) zf ,(g+; Q) 5(g++) T^{k^; P_) 5(g_)] , (3) 

where F^{Q'^) is the pion form factor and P± = P ± Q/2, q± = q ± P/2, qj^± = g_|_ ± Q/2 
and fc± = g ± (5/4. ^(g) is the dressed quark propagator, TT,{k; P) is the pion BS amphtude 
corresponding to relative qq momentum k (we choose equal partitioning) and pion momen- 
tum P, and Tr[. . .] denotes the trace over color, flavor and spin indices. The quark-photon 
vertex Ty{q\Q) at sufficiently large spacelike Q^ becomes Ty{q\Q) —>■ Q'j^, where Q is the 
quark charge operator. 

The quark-photon vertex satisfies the Ward-Takahashi identity (WTI) 

iQ,f^ip;Q)=Q{S-\p + Q/2)-S-\p-Q/2)) , (4) 

as a result of gauge invariance. At Q = the vertex is completely specified by the differential 
Ward identity 

^f,{p■,0) = Q^S{p)-'. (5) 

This reduces Eq. (^ to Ft^{Q'^ = 0) = 1 if Ftt is properly normalized^. This clearly shows that 
the bare vertex 7^ is generally a bad approximation due to the momentum-dependence of 
the quark self-energy: only with bare quark propagators, does a bare vertex satisfy Eqs. (^ 
and (^). Use of a bare vertex, in combination with dressed propagators, in Eq. (H), violates 
charge conservation and leads to P7r(0) 7^ 1. 

The WTI determines the longitudinal part of the vertex completely in terms of the 
(inverse) quark propagator. However, the transverse part is largely unconstrained by sym- 
metries. A common Ansatz for the dressed vertex is that due to Ball and Chiu |T^ which 
was developed in the context of QED investigations and is a representation in terms of the 
quark propagator functions A and B defined by 

S{pY^ = ti>A{p^) + B{p'). (6) 

With ff^ij); Q) = Q Vf^ij); Q) the BC Ansatz is 

rf (P;0) = H (-4(P.) + A(p_)) + 2fe /'^;)-^i''-) -2.p /"';)-g(^-' , (7) 

^ ^ P+ P- P+ P- 

where p± = p± Q/2. This satisfies the constraints from the WTI Eq. (^ and the Ward 
identity Eq. (^), transforms under CPT as a vector vertex should, and has the correct 
perturbative limit 7^ in the ultraviolet. The longitudinal part of T^'" is exact; the transverse 
part is exact only at Q = and in the UV limit. Curtis and Pennington have explored the 
additional Dirac structures that are possible for the transverse part of the vertex and have 



suggested an improved Ansatz [|1^] based on multiplicative renormalizibility in QED 



^This equivalence holds if the kernel of the pion BSE is independent of the pion momentum. For 
the ladder truncation of the kernel, which we consider here, this is the case. 



rr(p;Q) = rf(p;Q) + 

with M(p) = i?(p)/A(p). Both Ansatze Eq. (0) and Eq. (§) for the gg7 vertex satisfy all 
symmetry constraints, and use of them in Eq. (^) leads to -F,r(0) = 1, but neither of them 
contain the timelike vector meson poles of the exact vertex. 
The pion charge radius 



2 . QF^iQ^) 

ri = -Q 



dQ^ 



(9) 



will receive two types of contribution within the impulse approximation in Eq. (j^): 1) 
contributions from T^iQ = 0) coupled with the Q^-slope produced by the quark propagators 
and r^, and 2) contributions proportional to the Q^-slope of T^{Q). For the first type, use 
of the BC Ansatz is sufficient since the BC Ansatz is exact at Q = 0. However for the second 
type of contribution, (9r^(p; Q)/dQ^ is not constrained by symmetries. 

Previous studies of three-point quark loops such as Eq. (|^) for form factors have, for 
practical reasons, utilized parameterized representations of the DSE solutions for quark 
propagators and the BS amplitudes, in conjunction with the BC Ansatz for the quark- 
photon vertex UgHIOJ- The parameters are fitted to give a good description of pion and 



related chiral observables such as /jr, m.,^^ r.„-, and (qq) with a view towards parameter-free 
studies of other mesons and observables [Q. This procedure can produce values of r^ in 



the range 20% to 30% below the experimental value [[7|,p!4|, which leaves some room for 
additional contributions, such as those coming from vr vr rescattering [|l^ and from the p- 



resonance mechanism. However the estimated 50% of r^ being due to the p-resonance 
appears to be incompatible with these parametrizations. 

Here, we reconcile the VMD picture with the QCD picture of a photon coupled to 
distributed quarks in a qq bound state by calculating all ingredients needed for F^^IQ"^) via 
Eq. (^ from their dynamical equations: we solve the quark DSE, the pion BSE, and the 
inhomogeneous vertex BSE in a self-consistent way, using the same model for the effective 
quark- ant iquark coupling. Both resonance and non-resonance contributions to the vertex 
are dynamically generated through the BSE. By comparison with the BC Ansatz we can 
identify resonance contributions; however, one must keep in mind that such an identification 
is necessarily model-dependent. 

III. THE BETHE SALPETER SOLUTION FOR THE VECTOR VERTEX 

The quark-photon vertex satisfies the renormalized inhomogeneous BSE 

T,{p; Q) = Z2Ql, + j -^ K{p, q; Q)Siq + vQ)r,{q; Q)Siq - vQ) , (10) 

where t] + f] = 1 describes the momentum sharing between the two quarks. The kernel K 
operates in the direct product space of color, fiavor and Dirac spin for the quark and an- 
tiquark and is the renormalized, amputated qq scattering kernel that is irreducible with 



respect to a pair of qq lines. The notation j denotes a translationally-invariant regulariza- 
tion of the integral, with A the regularization mass-scale. At the end of all calculations the 
regularization is removed by taking the limit A — >^ oo. 

The renormalization constant Z2 and the renormalized dressed quark propagator S follow 
from the quark DSE 

^ 9'D,,{p - q)-^^^S{q)Tl{q,p) , (11) 

where D^^{k) is the renormalized dressed-gluon propagator and r"(g;p) is the renormalized 
dressed- quark-gluon vertex. The solution of Eq. (|l^) is renormalized according to 

S{p)-' ^ ^=t^ + mift), (12) 

at a sufficiently large spacelike /i^, with m{fi) the renormalized quark mass at the scale fi. 
In Eq. (|TT]), S, r^ and m{fi) depend on the quark flavor, although we have not indicated 
this explicitly. The renormalization constants depend on the renormalization point and the 
regularization mass-scale, but not on flavor: in our analysis we employ a flavor-independent 
renormalization scheme. 



A. Bound state contributions 

Solutions of the homogeneous version of Eq. ( p!OD at discrete timelike momenta Q^ define 
vector meson bound states with masses m^ = —Q^. It follows that T^{p]Q) has poles at 



those locations. The corresponding resonant form can be obtained by observing that Eq. (10 
has the equivalent form 



f,ip; Q) = Z2Qj, + Z2 j -^ M(p, q- Q)S{q + r,Q) Qj, S{q - rjQ) , (13) 

with M being the qq scattering amplitude given in schematic form by M = i^ + KSMS. 
In the vicinity of the timelike points Q^ = —m\ where the homogeneous BSE has solutions, 
M has poles with residues that define the physical meson BS amplitudes. In particular, for 

two flavors, 

Mha-c) rg(p;g)rg(g;-g) r-(p;g)r-(g; -Q) 

with the BS amplitudes properly normalized. In the ladder truncation that we will be 
concerned with in practice, the normalization condition reduces to (n = p^,uj) 

(15) 



^P^ = ^ll ^Tr [f ^(g; -K) S(g + r^P) T:{q- K) S{q - r^P)] 



P^=K'2= 



where the factor 1/3 appears because the three transverse directions are summed. Using 
the fact that massive vector mesons are transverse, the resonant form of the quark-photon 
vertex near the vector meson poles can be written as 

6 






where the couphng constants Qn {n = p , u) for p — 7 and a; — 7 mixing are 



mi Z2 r dl 



Qn S J {2tT 



Tr 



K{q;-Q)Siq + vQ)QluSiq-r^Q) . (17) 



Away from the pole, the separation into a resonant and non-resonant part is not unique; the 
solution of the BSE for the complete vertex contains both aspects in a consistent way. 

At the level of the ladder approximation, which is commonly used in practical calcula- 
tions, and using isospin symmetry, we have rUp = m^^ = my, and the flavor structure of the 
BS amplitudes is T^b; Q) = (ra/v^)!^^; Q), and r^(p;Q) = (1/72)1^^; <5)- The flavor 
trace in Eq. (|T7D gives g^^ = 3gp = 3gv- The factorization T^ = QT^ allows Eq. (|16D for the 
vector pole structure to simplify to the flavor-independent form 

where fyiTT'V = v^it^v / dv ■ Our numerical study is carried out for the flavor-independent 
vector vertex Tp{p; Q). 



B. Ladder rainbo^v truncation 

We use a ladder truncation for the BSE 

(\a \ru /■,„ \ts 

Y7mJ ®(y7.J , (19) 

where D^^l^ik) is the perturbative gluon propagator in Landau gauge. The resulting BSE is 
consistent with a rainbow truncation r"(g,p) — >■ 7,^A"/2 for the quark DSE, Eq. ([Tl|) , in the 
sense that the combination produces vector and axial- vector vertices satisfying the respective 
WTIs. In the axial case, this ensures that in the chiral limit the ground state pseudoscalar 
mesons are massless even though the quark mass functions are strongly enhanced in the 
infrared flJT6[|. In the vector case, this ensures electromagnetic current conservation. 



The model is completely specified once a form is chosen for the "effective coupling" Q{k'^). 
We employ the Ansatz @J^ 

^(P), (20) 



g{k'') _ 47r2 ^, 2^_fc2/,2 , ^^ 7mvr 



-DA;2e-^ /'^ + 47r 



iln 



r+(l + P/A2^CD)' 



with T{k'^) = [1 - exp(-fcV[4^t])]/^^ and 7^ = 12/(33 - 2A^/). This Ansatz preserves 
the one-loop renormalization group behavior of QCD for solutions of the quark DSE. In 
particular, it produces the correct one-loop QCD anomalous dimension of the quark mass 
function M{p^) for both the chiral limit and explicit chirally broken case P,p!7[|. The first 
term of Eq. (|20|) implements the strong infrared enhancement in the region fc^ = — 1 GeV^ 



which is a phenoinenological requirement for sufficient dynamical chiral symmetry breaking 



to produce an acceptable strength for the quark condensate |18[. We use m^ = 0.5 GeV, 
r = e — 1, Nf = 4, Aq^j-, = 0.234 GeV, and a renormalization point /x = 19 GeV, which is 
sufficiently perturbative to allow the one-loop asymptotic behavior of the quark propaga- 
tor to be used as a check 0,^. The remaining parameters are fixed to cj = 0.4 GeV and 
D = 0.93 GeV^ to give a good description of niT^/K and /,r- The subsequent values for fx 
and the masses and decay constants of the vector mesons p, 0, K* are very well described [^ . 

C. Numerical solution 

The general form of r^(g; Q) can be decomposed into twelve independent Lorentz covari- 
ants, made from the three vectors 7^, the relative momentum g^, and the photon momentum 
Qfj^, each multiplied by one of the four independent matrices 1, ^, Q, and crfj,uq^Qu- Four of 
the covariants represent the longitudinal components which are completely specified by the 
WTI and can be taken as the longitudinal projection F^ of the BC Ansatz, Eq. (|^). The 
solution of the BSE for the transverse vertex can be expanded in eight covariants T*(g; Q). 
Thus the total vertex is decomposed as 

8 

F^(g; Q) = Fj(g; Q) + J2 ^'M. Q) m\ Q " Q; Q') , (21) 

with the invariant amplitudes Fi being Lorentz scalar functions. The choice for the covariants 
T^(g; Q) to be used as a basis is constrained by the required properties under Lorentz and 
CPT transformations, but is not unique. The BSE, Eq. ([10|) , must be projected onto the 
covariant basis to produce a coupled set of eight linear equations for the invariant amplitudes 
Fi to be cast in matrix form. This requires a procedure to project out a single amplitude 
from the general form Eq. (^Tj). It is therefore helpful if the chosen covariants satisfy a 
Dirac-trace orthogonality property. The following set of orthogonal covariants is used here0 

T'^{q;Q)=l^,, (22) 

T,^(g;Q) = (gJf-|7J(gT)/g^ (23) 

T',iq;Q) = q^,Qq-Q/iq'Q'), (24) 

T^iq^Q) = - KlQ, ^] +2gJ(^) /2g, (25) 

T',{q;Q)=tqj:/q, (26) 

T',{q;Q) = ^bj:,f]g-Q/<i', (27) 

Tliq; Q) = ^ bl m (1 - cos^ e) - 2 T^«(g; Q) , (28) 

T'^{q■,Q)=^qlfQ/q\ (29) 



^Since the domain of interest here includes both timelike and spacelike Q^, we avoid factors of 
Q = VQ^- Therefore we do not require the covariants to be normalized, as we did in Ref. [^. We 
also include explicit factors of q ■ Q mT^ and Tq so that now every Fi is even in q ■ Q. 



where V'^ is the component of V transverse to Q, that is V^ = V^ — Q^{Q ■ V)/Q'^ and 
q ■ Q = qQ cos6. Each amphtude can be projected onto a basis of even-order Chebyshev 
polynomials in cos6' to produce "Fj(g^,(5^) 0. It is found that terms beyond the zeroth 
order, n > 0, are insignificant and are not used in the subsequent study of the form factor. 

For a range of small timelike and spacelike Q^, the ratio of the BSE solution to the 
BC Ansatz for the dominant amplitude °Fi is displayed in Fig. |^ as function of g^. The 
deviation from unity indicates differences between the BC Ansatz and our numerical solu- 
tion of the ladder BSE. It is immediately obvious that for timelike photon momenta and 
g^ < 3 GeV^ there is a significant enhancement of the BSE solution compared to the BC 
Ansatz. This strong increase in the vertex BSE solution is due to the vector meson pole 
at Q^ ^ —0.55 GeV^, which only affects the vertex at small g^; at large g^ both the BSE 
solution and the BC Ansatz approach the bare vertex for any value of Q^. 

Charge conservation is evident through the ratio being one at Q = 0: the ladder vertex 
BSE solution satisfies the WTI. The fact that our numerical solution indeed is to within 
1% equal to the BC Ansatz, not only for °Fi, but also for the other amplitudes, indicates 
the accuracy of our numerical methods. At small spacelike momenta the BSE solution 
differs considerably from the BC Ansatz. The difference increases with Q^ in the range 
< Q^ < 2 GeV^; it is only at asymptotic values of Q^ that they become equal again. 

A different perspective is shown by Fig. ^ which compares the Q^-dependence of ^Fi at 
g^ = 0. Again, the agreement at Q^ = is dictated by the Ward identity; however, the slope 
of "^-^1(0, Q"^) aX Q^ = as obtained from the vertex BSE solution differs significantly from 
that produced by the BC Ansatz. This can be of importance for the charge radii of the pion 
and other hadrons. In this figure one can clearly see the presence of the vector pole in the 
gg7 vertex at Q^ ~ —0.55 GeV^. Also in this figure we plot the amplitude "Fi as obtained 
from the BC Ansatz with the quark propagators of Ref. 0], which were parameterized 
rather than obtained as the solution of a DSE. The parameters were fitted to f^^, m^r, r^^, 
and (gg), using the BC Ansatz in the calculation of r^^. The Q^-dependence of the gg7 
vertex necessary in order to reproduce a reasonable value for r^^ was thus incorporated via 
the quark propagator functions. We see that this phenomenological vertex amplitude has 
indeed a much stronger Q^-dependence than is obtained by use of the present DSE solution 
for the quark propagator in the BC Ansatz. This phenomenological procedure allows one to 
get reasonable electromagnetic form factors at low Q^ without explicitly taking resonances 
into account. 

We have investigated the extent to which the vertex BSE solution can be represented by 
the interpolating form 

r,(g; Q) ^ rj^(g; Q) -J^ T;(g; Q) FY {q') ^Jq,^_^^.^^ , (30) 

where Fy are the leading Chebyshev moments of the mass-shell vector meson BS amplitudes 
produced by the present model |^. The second term Eq. (|30D correctly describes the vertex 
near the vector meson pole at Q^ = —my] there is no width generated for the vector meson 
in the present ladder truncation of the BSE. One would have to add the tttt channel to the 
ladder BSE kernel to produce a p width; for the vertex, this would generate an imaginary 
part beyond the threshold for pion production, Q^ < — 4m^. 



For small spacelike Q^, Eq. ( pOD represents an extrapolation consistent with: a) 
r^^(g; Q = 0) is exact due to gauge invariance and b) the p — 7 coupling should not generate 
a photon mass. Over the hmited domain —my < Q'^ < 0.2 GeV^, Eq. (|30| ) produces a very 
good representation of the vertex amplitudes °Fi, "^3 and ^F^ and a quite reasonable repre- 
sentation for the remainder, except for ^^4 which is clearly the worst case, and the hardest 
to represent in a resonance formula. In Figs. |^, ^ and P, this approximate form is compared 
to the vertex BSE solution for the three dominant amplitudes °Fi, ^^4 and ^F^. 

In this sense we shall refer to the BC Ansatz as missing a contribution near the 
p pole, and to Eq. (^) as a resonant-improved BC Ansatz applicable in the region 
—0.4 GeV^ < Q^ < 0.2 GeV^ only. At asymptotic Q^, this form is obviously not adequate: 
it does not go to the bare vertex. Close to the p pole, neither our vertex BSE solution 
nor the interpolating form Eq. (^) are physically realistic, because neither one takes into 
account the width Tp generated by the open vrvr decay channel. Since Vp/rrip ~ 0.2, the 
width contributes little for Q^ > —0.4 GeV^, and the p contribution can be approximated 
by the simple pole term in Eq. (^) for applications such as the pion form factor calculation. 

IV. RESULTS FOR THE PION FORM FACTOR 

In addition to the solution of the BSE for the gg7 vertex and the quark propagator S as 
the solution of its DSE, we also need the pion BS amplitude TT^lq; P). This BS amplitude is 
the solution of a homogeneous BSE, and for pseudoscalar bound states it has the form Q] 

r,(g; P) = 75 [lE^iq; P) + f F^{q-P) + iq- P G^{q; P) + a^, q^P, H^{q; P)] , (31) 

with the invariant amplitudes E,^, F^, Gj^ and H,^ being Lorentz scalar functions of g^ and 
q ■ P = qP cos 6, with P^ = —rri^ the fixed, on-shell, pion momentum. Each amplitude can 
be projected onto a basis of even-order Chebyshev polynomials in cos 6 and in Ref. [Q it was 
found that terms beyond the zeroth order are insignificant. With these ingredients we can 
now calculate the pion electromagnetic form factor in impulse approximation, Eq. (|^). 

In Fig. 1^ we show our results for the pion form factor for a range of timelike and space- 
like Q^ using the BSE solution for the qq'y vertex, and compare our results with various 
other treatments of this vertex. Clearly a bare quark-photon vertex is incorrect: current 
conservation, which ensures -F7r(0) = 1, is violated. Use of the BC Ansatz conserves the 
current because it satisfies the vector Ward identity. However the resulting form factor 
misses the data completely; the Q^-dependence is clearly too small. On the other hand, 
the BSE solution for the vertex gives an excellent description of the low Q^ data, both in 
the spacelike and timelike region, without fine tuning the model parameters: the parameters 
are completely fixed by m^ and /^ in Ref. P). The use of the Curtis-Pennington Ansatz, 
Eq. (^), gives results for the form factor essentially the same as those obtained with the BC 
Ansatz. 

Experimentally, the form factor shows a resonance peak at Q^ = —m'l = —0.593 GeV^. 
Our calculated Ft^^Q"^) using the BSE solution diverges as l/iQ"^ + my) as one approaches 
Q'^ = —my ^ —0.55 GeV^, where the homogeneous BSE admits a bound state solution: in 
the ladder truncation for the BSE one does not get a width for the p meson. We expect 
that if we incorporate the open decay channel p — >■ tttt in the BSE kernel to generate a 

10 



width Tp for the p, we will have better agreement with the data close to the resonance 
peak. In addition, some of the present difference between our calculation and the data in 
the resonance region is due to the fact that our numerical value for rrip is about 4% below 
the physical p mass. 

Four of the eight amplitudes, F3, Fg, F7, and Fg, contribute less than 1% to -FV(Q^) 
on the Q^-range considered. Although the amplitude F3 contributes very little to the form 
factor, this amplitude is needed for strict charge conservation. Therefore we can truncate the 
BSE solution to the dominant five covariants Fi through F5 without significantly changing 
the description of F^(Q^). In Fig. |^, the result using the five dominant amplitudes is almost 
indistinguishable from the results obtained with all eight amplitudes. 

The relative contributions to Ft^{Q'^) from the amplitudes Fi, F2, F4, and F5 of the vertex 
BSE solution are displayed in Fig. ||. This shows that the canonical amplitude Fi, together 
with the canonical (pseudoscalar) amplitude F^ of the pion, generates the bulk of the form 
factor in the infrared region. The pion pseudovector amplitudes F^ and G^^ give a negative 
contribution of about 25% near Q^ = 0, but they are known to dominate at large Q^: the 
asymptotic behavior of the form factor, F^^IQ^) ~ 1/Q^, is governed by the pseudovector 
amplitudes of the pion |Q. The vertex amplitudes F2, F4, and F5 give negative contributions 
to Ftt of the order of 10% at small Q^. Of particular note is the contribution from amplitude 
F4. Due to the Ward identity, this vertex contribution, and hence its contribution to F^^^Q"^), 
must vanish at Q^ = 0. The latter is indeed evident in Fig. |; it is also evident there that 
this amplitude nevertheless is one of the important contributors to the slope and hence r^r; 
its contribution to the form factor grows roughly linear with Q^ in the range we studied. 
This could have significant consequences for the form factor at intermediate energies. At 
asymptotic values of Q^ the gg7 vertex approaches the bare vertex, and thus all functions 
Fi except for Fi will vanish. The asymptotic behavior of the form factor is therefore not 
influenced by F4 nor any other subdominant amplitude of the gg7 vertex. 

In terms of the pion charge radius r^r, these results are summarized in Table |I[ It is 
seen that, compared to the empirical value r^ = 0.44 fm , the BSE solution for the vertex 
generates an excellent value: r^ = 0.46 fm^ using all eight transverse covariants or using the 
dominant five, Fi through F5. The BC Ansatz produces a value for r^ that is less than half 
this value; this is also evident in Fig. |^. The main reason for this is that the gg7 vertex 
at small but non-zero Q^ is poorly represented by the BC Ansatz, as shown, for example, 
in Fig. ^. The Curtis-Pennington Ansatz gives almost the same charge radius as the BC 
Ansatz. 



A. Resonance contribution to r^^ 

In the previous section we found that the vertex BSE solution could be simulated by a 
resonant-improved BC Ansatz, Eq. (^). The use of such a form for the gg7 vertex yields, 
via the impulse approximation Eq. (H), a charge form factor that can be expressed as 



F^) ^ F^'^iQ') - ^P- ;^;;-^^'\^' , (32) 



5fp(Q2 + m 



where Ff^((5^) is the result from the BC Ansatz. The combination gp^^Fy-Kn represents 
the TTTT coupling to the vector qq correlation T'lF^ evident in the resonant term of Eq. (pOD. 
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The behavior of this latter term guarantees that at the vector meson pole, the second term 
of Eq. (|32D reproduces the standard behavior and FvTTn{—^p) = 1- For other Q^ values 
Fv-KniQ'^) does not correspond to a physical process: off-shell mesons are by definition not 
physical. However, the departure of Fv-KniQ^) from unity is a rough measure of the difference 
in TTTT coupling experienced by the effective vector qq correlation away from the p mass-shell 
compared to the physical pTrvr coupling. 

The charge form factor resulting from two versions of the resonant-improved BC Ansatz 
are shown in Fig. p. One version is based on the five dominant invariant amplitudes ^Fy{q'^) 
of the p meson, the other uses only the single dominant amplitude ^F^iq^)'-, physical nor- 
malization is imposed in each case. The employed values of m^ and fp in each case are those 
obtained consistently from the homogeneous BSE in the present model [^ , see Table |I|. It is 
seen from Fig. ^ that both forms of the resonant-improved BC Ansatz simulate the behavior 
of -F^(Q^) produced by the vertex BSE solution quite well. 

The obtained values of r^ are shown in Table | for the approximations discussed above. 
Also shown there is the charge radius produced by using the two dominant amplitudes Fi 
and -F5 for the resonant addition to the BC Ansatz. Compared to r^BC ~ O-^^ fm , all 
versions of a resonant-improved BC Ansatz provide results roughly a factor of two higher 
and in significantly better agreement with the result from the vertex BSE solution and the 
experimental data. When the five dominant p amplitudes are used, the resonant-improved 
BC Ansatz reproduces the vertex BSE result for r^ to within 5%. The error increases to 
10% if only the single dominant vector meson covariant is used for the resonant term of the 
vertex; part of the reason for this is due to the increase in the vector meson mass : with 
Fi only, rn? = 0.766 GeV^, which moves the pole 30% further into the timelike region. 

Use of Eq. (^) allows r^ to be characterized as a resonant addition to the result of the 
BC Ansatz. If we write the difference between r^ and the BC contribution r^ g^ as 

and compare this with Eq. @, we see that -FV7r7r(0) characterizes the necessary weakening 
of the VMD mechanism for r^ to account for the distributed qq substructure. The values 
of FvnniO) defined by Eq. (|33D , using the experimental values for gpmr, rup, and Qp, are 
also shown in Table |. The difference between the r^ result from the vertex BSE solution 
and r1 g^ represents about 60% of the VMD value for the charge radius, r^ ^y^-^ = 0-48 fm^. 
The result from the best resonant-improved BC Ansatz (using the vector meson amplitudes 
Fi — -F5) shows that the resonant term is contributing to r^ at the level of 50% of the VMD 
value. This decreases somewhat when the description of the vector meson is simplified. One 
can attribute this weakening of the VMD mechanism at Q^ = to the fact that the photon 
couples to a distributed, interacting qq correlation and a significant part of this is already 
accounted for by the BC Ansatz for the vertex. With the remainder viewed as due to the p 
resonance, its effect is overestimated if, for example, the coupling to vrvr is described by the 
physical mass-shell value gpT,T,. The effective reduction of the gp.„T, is evidently about 50% in 
the present approach; this supports the previous cruder estimate [0]. 

In a more phenomenological approach, using parameterized quark propagators fitted to 
pion observables including r^^ as was done in Ref . [0 , one obtains a much better value for r^ 
using a BC Ansatz. In such an approach the Q^- dependence of the vertex is parameterized 

12 



via the quark propagator, and can thus give a fair representation of the pion form factor. 
The charge radius obtained in this way, r^ = 0.55 fm, is much closer to the experimental 
value than the one in the present work using a BC Ansatz. This is also an indication of 
the model-dependence and arbitrariness of separating the contributions into resonant and 
non-resonant terms. 

The utility of the BC Ansatz for studies of electromagnetic coupling to hadrons is that it 
is completely specified in terms of the quark self-energy amplitudes. For similar reasons we 
seek to summarize the main features of the BSE solution for the quark-photon vertex also in 
terms of the quark self-energy amplitudes and a phenomenological p-meson BS amplitude, 
based on the resonance formula Eq. (|30|). For the p BS amplitude we use a simple Ansatz 
for the dominant amplitude only, FY{q^) -^ ^p/{^ + Q^/^'^)-, properly normalized and with 
uj = 0.66 GeV to give a decay constant /p = 201 MeV, similar to the parameterization used 



in Ref. |]14|. In order to construct an Ansatz which can be used both at small and at large 



Q^, we consider the following phenomenological form for Q^ > —m?p 

i^W.^J i;. W,^j % l + qy^^rn.iQ^ + ml)^ ' ^"^^^ 

Since the correct UV asymptotic limit is provided by rl^*-^, the parameter a > provides for 
a Q^ suppression of the resonant term. We fit a to r^, using the experimental values for rup 
and fp, and find good agreement with a = 0.03 GeV^^. The form factor F^i^Q"^) produced 
by this phenomenological vertex is shown in Fig. ^, and is seen to provide an excellent 
fit to the vertex BSE solution, and hence according to Fig. ^ also to the experimental 
data. Extensions to a more realistic p BS amplitude are straightforward to implement in 
phenomenological model calculations. 



V. SUMMARY AND OUTLOOK 

We have used a ladder truncation for the BSE for the quark-photon vertex, in conjunction 
with a ladder-rainbow truncation for the quark DSE. Both the vector WTI for the quark- 
photon vertex and the axial-vector WTI are preserved in this truncation. This ensures both 
current conservation and the existence of massless pseudoscalar mesons if chiral symmetry 
is broken dynamically: pions are Goldstone bosons [^ . The details of the model were fixed 
in previous work 0. It leads to dynamical chiral symmetry breaking and confinement; fur- 
thermore, at large momenta, our effective interaction reduces to the perturbative running 
coupling and thus preserves the one-loop renormalization group behavior of QCD and re- 
produces perturbative results in the ultraviolet region. The model gives a good description 
of the TT, p, K, K* and masses and decay constants 0. 

Our numerical solution of the vertex BSE shows clearly the vector meson pole in all eight 
transverse amplitudes. At the photon momentum (5 = 0, the solution agrees perfectly with 
the BC Ansatz, as required by the WTI and gauge invariance. Also at spacelike asymptotic 
momenta our transverse solution agrees with the BC Ansatz: both go to the bare vertex. 
At small but nonzero Q^ the BSE solution departs significantly from the BC Ansatz: the 
BSE solution has a stronger Q^- dependence. In the region —m? < Q^ < 0.2 GeV^, the 
quark-photon vertex can be described by a BC Ansatz plus a resonant term. However, there 
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is no unique decomposition of the vertex into resonant and non-resonant terms away from 
the pole, and the BSE solution for the vertex is the appropriate representation containing 
both aspects. 

Subsequently, we have calculated the pion charge form factor in impulse approximation. 
The results using the BSE solution for the gg7 vertex are in excellent agreement with the 
data, see Fig. |^, and produce r^ = 0.46 fm , without fine tuning the model parameters: the 
parameters are completely fixed in Ref. |^. Use of the BC Ansatz generally leads to a charge 
radius which is too low: in the present model, the resulting r^ = 0.18 fm^ is less than half 
the experimental value, while the remainder can be described by a p resonant term. This 
indicates that as much as half of r^ can be attributed to a reasonable extrapolation of the p 
resonance mechanism. On the other hand, the strict VMD picture is too simple; about half 
of r^ arises from the non-resonant photon coupling to the quark substructure of the pion. 
One should however keep in mind that such a separation in resonant and non-resonant 
contributions is ambiguous and model-dependent. 

The form factor Fj^{Q'^) exhibits a resonance peak at timelike momenta Q^ near — m^, 
and our calculated Ft^^Q"^) does indeed show such a peak. Since in ladder truncation one 
does not generate a width for the p meson, we overshoot the data close to the p-pole. We 
expect that if we include the tttt mechanism for the p width in our formalism, we will have 
better agreement with the data in the timelike region. A detailed comparison close to the 
p-pole will be postponed until we have included this effect in our calculations. 

Pion loops will not only generate a width for the p, they will also give a direct contribution 
to the pion form factor and generate a nonzero imaginary part for momenta Q^ < —Am"^. 
Estimates are that the pion charge radius increases by 10% to 15% |15[ due to pion loops. 



which would lead to a charge radius which is too large in the present model. However, the 
pion loops will also affect the calculated meson masses: estimates for the shift in the p mass 
vary between 2% and 10% ||2^, and the parameters in the effective quark- ant iquark coupling 



will have to be readjusted to maintain agreement with the experimental data. We hope to 
address these question in future work. 

With our BSE solution for the qq'y vertex we can now investigate other hadron form 
factors. First results indicate p3[ that also the 7*7r7 transition form factor at low Q^ is 
reasonably well described with the present vertex BSE solution. In the near future we 
will also apply this method to the kaon form factor and to electromagnetic decays such as 
p ^ TT'y. Another interesting application is the nucleon form factor, using models for the 
nucleon such as those being developed in Refs. [0,B3]. The phenomenological resonance 



addition to the BC Ansatz, Eq. (|34D, can be used as guidance for phenomenological models 
to facilitate these and other hadron form factor calculations, without having to solve the 
vertex BSE numerically. 
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FIGURES 




FIG. 1. The impulse approximation for the pion charge form factor. 
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FIG. 2. The ratio of the BSE solution over the BC Ansatz for the dominant amplitude ^Fi, 
associated with 7J, of the quark-photon vertex as function of q"^ for the indicated values of the 
photon momentum Q^. The quark momenta are q it Q/2. 
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FIG. 3. The Q^-dependence of the BSE solution for the dominant amphtude ^Fi of the ^g7 
vertex at g^ = (soHd hne). The dashed hne with circles is the BC Ansatz for this amplitude in 
the present model; the dotted line with diamonds is the BC Ansatz with phenomenological quark 
propagator functions A and B parameterized and fitted to pion observables, including r^^, according 
to Ref. 0. 
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FIG. 4. The g^-dependence of the BSE solution for the dominant amplitude ^Fi for a range 
of timelike and spacelike values of the photon momentum Q^. The circles correspond to the 
resonant-improved BC Ansatz as described in the text. 
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TABLES 



photon vertex 


rn 


rl 


2 2 

'''n ~ ^7r,BC 


Fv.AO) 


Expt. H 


0.663 ± 0.0006 


0.44 ± 0.0012 






DSE Fi - Fg 


0.678 


0.460 


0.278 


0.58 


DSE Fi - F5 


0.677 


0.459 


0.277 


0.58 


BC Ansatz (BCA) 


0.426 


0.182 






EGA + res. Fi - F5 


0.640 


0.409 


0.227 


0.47 


BCA + res. (F1&F5) 


0.625 


0.390 


0.208 


0.44 


BCA + res. Fi 


0.605 


0.366 


0.184 


0.38 


BCA + phenom.res. Fi 


0.68 


0.46 


0.28 


0.58 


BCA of Ref. § 


0.55 


0.30 







TABLE I. Results for the pion charge radius according to various treatments of the 
quark-photon vertex. FvmriO) roughly characterizes the reduction in the vrvr coupling strength 
for the vector qq correlation in the vector vertex compared to the on-shell coupling to the p, gpmr- 





mp 


fp 


Expt. 


0.770 


0.216 


Fi-Fg 


0.742 


0.207 


F1-F5 


0.730 


0.201 


F1&1F5 


0.765 


0.203 


Fi 


0.875 


0.20 



TABLE IL The vector meson mass rup and decay constant fp in GeV corresponding to various 
truncations of the eight possible invariant amplitudes of the p BS amplitude [0. 
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